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Abstract 

O Using the pole approach we determine the mass and width of the 

/o(980), in particular we analyze the possibility that two nearby poles 
D ! are associated to it. We restrict our analysis to a neighborhood of 

the resonance, using tttt data for the phase shift and inelasticity, and 
the invariant mass spectrum of the J/ip — * (ftirir,4>KK decays. The 
formalism we use is based on unitarity and a generalized version of the 
Breit-Wigner parameterization. We find that a single pole describes 
the /o(980), the precise position depending upon the tttt data used 
(see Eq. (fTU|) of the main text). As a byproduct, values for the 
and gf KK coupling constants are obtained. 



1 Introduction 



According to the Review of Particles Physics (RPP) pQ, the mass of the 
/o(980) scalar resonance is 980 ± 10 MeV whereas the width ranges from 
40 to 100 MeV. The reasons for such an uncertainty in the width are the 
great amount and variety of experimental data and the different approaches 
used to extract the intrinsic properties of the resonance. To these points 
we could also add the lack of a precise definition of what is meant by mass 
and width, although there seems to be consensus in using the pole approach, 
where the mass and width of the resonance are found from the position of the 
nearest pole in the T-matrix (or equivalently, the S'-matrix). However, even 
if this approach is adopted the final results differ on the number, location and 
physical interpretation of the poles. This is because the pole approach is not 
enough to completely fix the framework needed to perform the resonance 
analysis; in fact there are different formalisms that made use of it. An 
example is field theory, where a finite imaginary part of the propagator arises 
after Dyson summation of the one-particle-irreducible diagrams contributing 
to the two-point function. Unitarity also implies a general complex structure 
of the T-matrix in terms of which the pole approach becomes relevant to 
define the mass and width of the resonance. However, the general solution 
to the unitarity constraint has no implications regarding the number and/or 
locations of the poles. Thus, most analyses using the pole approach must 
involve further assumptions. 

Far from physical thresholds, the identification of the mass and width of 
a resonance in terms of the nearest pole in the T-matrix is not ambiguous. 
However, when the resonance lies in the vicinity of a threshold this identifi- 
cation is not so obvious and more than one single pole can be required for 
a correct description of the resonance (see for example Ref. This could 
be the case for the /o(980) whose mass is very close to the KK threshold. 
In Ref. jH], Morgan and Pennington (MP) use a formalism general enough 
to avoid any assumption about the number of poles associated to a reso- 
nance. For the particular case of the /o(980), their exhaustive analysis leads 
to the conclusion that the / (980) is most probably a Breit-Wigner-like res- 
onance — with a narrow width T ~ 52 MeV — which can be described in 
terms of two nearby poles (in the second and third sheets). Moreover, the 
precision data coming from J/ip — > <fi(MM) decays, where (MM) stands for 
7T7T or KK, play an essential role as a crucial check in favour of the two-pole 
description of the /o(980) and disfavour the cases with one or three poles. 
These results indicate that the description of the /o(980) using a Breit-Wigner 
parametrization seems to be appropriate and should give results similar to 
those in Refs. |3J E]. However, an existing analysis using a Breit-Wigner 
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parametrization performed by Zou and Bugg (ZB) concludes that the 
/o(980) is most likely a resonance with a large decay width (~ 400 MeV) 
and a narrow peak width (~ 47 MeV). Later on, MP jH] and ZB [6. have 
both confirmed their former results, leaving the agreement between the two 
approaches as an open question. 

The pole approach formalism has been successfully applied to hadronic 
resonances such as the p(770) [3 El El and the A(1232) [TU] . An important 
advantage of this formalism is that it yields process and background indepen- 
dent results. This background independence is only valid if other resonances 
are not present within the kinematical region under consideration 1 . Thus, if 
one insists on this point, as we will, it is important to restrain the analysis 
to a neighborhood of the resonance under study. Due to the previous consid- 
eration, we exclude from our analysis of the /o(980) the central production 
data pp — > p{tttt,KK)p [H]-[IE], which covers a much wider energy range 
and whose phenomenological description requires not only tttt and KK scat- 
tering but also a production mechanism involving many parameters. We also 
exclude the experimental data on J/t/; -> uj(nn, KK) decays because the f 
signal is too weak fT7\ ITS] 2 and the inclusion of these data in the fit would 
require further parameters 3 . 

Our purpose in this work is to extract the mass and width of the /o(980) 
scalar resonance using the pole approach. As a byproduct, we also obtain 
values for the coupling constants g/ own and Qf KK- We suggest that, for nar- 
row resonances as the /o(980), the complete one-loop scalar propagator must 
be used in the pole equation. Thus, our analysis is based on a generalized 
Breit-Wigner description of a scalar resonance coupled to two channels, not 
only satisfying unitarity but also including loop effects. As far as the pole 
is concerned, we pay special attention to the possibility of describing the 
/o(980) in terms of more than one pole, a phenomenon which is known to 
occur when the mass of the resonance is close to a threshold. Nevertheless, it 
is worth remarking that the need of two poles is not guaranteed, it strongly 
depends upon the precise value of the renormalized mass of the resonance 

lr The amplitude associated to a given resonance is not expected to describe the physics 
in a large kinematical region (compared to the width of the resonance) where additional 
resonances can exist. 

2 After submission of this manuscript the BES Collaboration has confirmed these results 
PI 

3 A phenomenological analysis of the J/tp — > VS decays assuming U(3) symmetry 
predicts that for ideally mixed scalar and vector states the values of the relevant coupling 
constants are: gj/^h = 9j/^a = 1 and g.j/^f = = 0. The small departure from 

ideal mixing of the <fr and the /o states explains, in a first approximation, the importance 
of the fo contribution to the J/ip — > 4>(tttt,KK) decays and, on the contrary, the minor 
role played in the J/ijj — > ui(tttt, KK) decays 
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and its coupling constants to the two channels. 

Values of the renormalized mass tur and the coupling constants gf mr and 
9 foKK are obtained from a fit to experimental data including the nir phase 
shift and inelasticity as well as the J/ip — » (f)7nr(KK) invariant mass spectra. 
We then look for poles in the four Riemann sheets associated to a resonance 
coupled to two channels. Our conclusion is that the / (980) can be described 
in terms of a single pole whose precise position depends upon the tttt data 
used (see Eq. (fT9j) of the main text for details). 



2 Formalism 

Before discussing the formalism used for the particular case of the / (980), 
it is convenient to briefly summarize two well-known definitions of mass and 
width of a given resonance, both widely used in the hadron physics literature 
(see Ref. JT0] and references therein). One definition, known as the conven- 
tional approach, is based on the behavior of the phase shift of the resonance 
as a function of the energy, while the other, known as the pole approach, is 
based on the pole position of the resonance, which as discussed in the intro- 
duction includes several approaches. We will not consider here the powerful 
formalism developed in Ref. [Sj since it goes beyond a Breit-Wigner-like de- 
scription of the resonance, to which our analysis is restricted. In this sense, it 
is worth noticing that these more powerful methods provide further support 
to the Breit-Wigner description of the /o(980). Thus, the analysis carried in 
this paper is more restricted in scope, although it turns out to be general 
enough for the /o(980) case. 

In the conventional approach, the mass and width of the resonance are 
defined in terms of the phase shift 5 as 4 



5{s = M\ 



90° 



1 

M~ s 



d5{s) 
ds 



-i 



(i; 



respectively. Since the phase shift is extracted from direct comparison with 
experimental data, the decay width defined in this way is usually called the 
visible or peak width. For an elastic Breit-Wigner (BW) resonance the phase 
shift is chosen as 5 

MsTs(s) 

tan5( S ) = --- w , (2) 



4 We use the subindexes 8 and p to denote the mass and width of the resonance in the 
conventional and pole approaches respectively. 

5 A s-dependent width ra(s) in Eq. (J2J) is mandatory when a background around the 
resonance is assumed [TO] . 
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which leads to the partial-wave amplitude 



= e 2lS - 1 MsTs(s) 
a 2i s-Mf + iMsTs(s) ' 1 } 

where s is the center-of-mass energy squared. 

In the pole approach (or ^-matrix approach), the resonance shows up as 
a pole in the amplitude 

a = — ^— + B , (4) 

S Sp 

where the two terms correspond to the resonant and background contribu- 
tions separated according to Refs. [2U 122] ■ Eq- (J1J) is understood as a power 
series expansion of the amplitude around s p , therefore, in order this descrip- 
tion to make sense, the background around the pole (which is fixed from the 
fit to experimental data) should be a smooth function of s affecting minimally 
the pole position. In this approach, the mass and width of the resonance are 
defined in terms of the pole position s p as 6 

Sp = m 2 p - irripTp . (5) 

The pole approach provides a definition for the parameters of an unstable 
particle which is process independent (independent of the process used to 
extract them) and also background independent (different parametrizations 
of the background will hardly modify the values obtained for the pole pa- 
rameters of the resonance) . 

In the remaining of this section we pursue the pole approach for the case 
of a resonance coupled to two channels, including furthermore the possibility 
of a strongly s-dependent width due to the opening of a second two-body 
threshold. Two ingredients are required in order to build the scattering 
amplitude to be used in our formalism: unitarity and the complete one-loop 
scalar propagator. 

Concerning the first of the ingredients, unitarity sets stringent constraints 
on the amplitudes needed for the description of a resonance coupled to several 
channels. The correct incorporation of these constraints into the S"-matrix is 
compulsory for an adequate analysis of experimental data. The analysis in 
the general case would require a model independent approach, as in Ref . |3] , 
to determine the number and location of poles associated to the resonance. 
However, previous analyses [3 E] have shown that the / (980) can be de- 
scribed in terms of a Breit-Wigner-like resonance with two poles associated 



6 The relationship between the mass and width parameters defined in the conventional 
and pole approaches can be found in Ref. [T77] . 
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to it. The Breit-Wigner parameterization is nevertheless a particular case 
of an amplitude fulfilling unitarity. Indeed, for a relativistic particle, the 
general solution to the unitarity constraint can be written as j22j 



_ e 2iS °& - 1 e *(*aW+*6(«)) ^/G^Gl 
ah ~ ab " F(s) + iG(s) ' (6) 

where 8 a b stands for the phase shifts describing the background in channels 



(a, b), P a ,b = yl -4m^/s with m aj b the masses involved in the two-body 
decays of the two channels, F(s) and G(s) are arbitrary real functions of 
s, and Gi(s) are positive functions with the property G(s) = J2i= a ,bGi(s). 
Identifying F(s) = s — mi and G(s) = m p T p , the amplitude © reduces in 
the one channel case to the amplitude (jlj) up to an overall normalization 
factor. 

The second of the ingredients mentioned above will be used in our frame- 
work in order to identify the functions F(s) and G(s) to the real and imag- 
inary parts of the complete one-loop propagator respectively. The previous 
identification has the advantage of incorporating automatically thresholds ef- 
fects (see below). This procedure requires the use of an effective field theory 
in order to calculate the full propagator of the / (980). In general, effective 
field theories are of limited use in the description of hadron physics where 
one expects the interactions to be strong. There are cases, however, where 
these theories can be used. The treatment of width effects, when the width to 
mass ratio (taken as an expansion parameter) is small, is an example where 
effective field theories can be useful for the description of narrow resonances, 
but not for broad ones. Notice in this respect that for the /o(980) Tf /mf ~ 
0.04-0.1 Pp. Concerning the final form of the scalar propagator, the use of 
a simple Breit-Wigner parametrization with constant width, which is appli- 
cable only to narrow resonances far from thresholds, is not enough due to 
the closeness of the KK threshold and the fo mass. Instead, one could use 
an energy dependent width, incorporating the kinematic dependences on the 
energy, but this approach amounts to include only the imaginary part of the 
self-energy. In our analysis, we prefer to use the fully corrected one-loop 
propagator, including both the real and imaginary parts of the self-energy, 
since this approach allows us to have a consistent description of the analytical 
properties, i.e. it provides a proper analytic continuation of the scattering 
amplitude below the KK threshold. 

After Dyson summation, the propagator of a scalar particle is 



p 2 — ttlq + lI(P 2 ) 
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where m is the bare or tree- level mass of the resonance and II (p 2 ) is the 
one-particle-irreducible (1PI) two-point function. In the on-shell scheme, a 
Taylor expansion of the real part of II (p 2 ) around the resonance mass allows 
to rewrite the scalar propagator as 

where the renormalized mass tur (the so-called on-shell mass) and the wave- 
function renormalization factor Z are defined as 

m 2 R = m 2 - ReU{m 2 R ) , 

Z- 1 = 1 + RelT(m|) , " ' 

with RelTQo 2 ) = dKeU(p 2 ) / dp 2 . By analogy with a Breit-Wigner resonance, 
the width is defined by 7 

T R = -^-ZlmU(m 2 R ) . (10) 

However, this on-shell definition of the resonance width is inadequate since 
it vanishes when a two-particle s-wave threshold is approached from below 
Due to the failure of the Taylor expansion of IT(p 2 ) around m 2 R , Eq. (fTUJ) 
has not the desired behaviour for a width properly defined. This is precisely 
the case under consideration since the KK threshold lies in the vicinity of 
the / (980) mass. 

On the contrary, the pole approach provides a consistent definition of 
the resonance width that behaves sensibly in the threshold region. In this 
approach, the Taylor expansion of II (p 2 ) is not performed and the scalar 
propagator (J7J) is written as 

A(p 2 ) = p2 _ m g t + n( p 2) 

~~ p 2 -m 2 R +ReU(p 2 )~RcU(m 2 R )+ilmU(p 2 ) ' 

where m 2 R = j«q — KeU(m 2 R ). Within the framework of the pole approach, 
the scattering amplitude (JHJ) describing a resonance coupled to two channels 
is obtained by identifying the functions F(s) and G(s) with the denominator 
of the complete one- loop scalar propagator in Eq. (fTTj) : F(s) = s — m 2 R + 
Rell(s) — ReII(m|j) and G(s) = ImII(s). This procedure leads to 

ab ~ 2i^K° ab S -m 2 fl +Ren(s)-Ren( m 2 J )+iImn( S ) 
~ 2iy/plp~ b ab 16tt s-m^+Ron(s)-Rcn(m^)+iImn( S ) ' 

7 This definition applies only to narrow resonances, <C to/j, where ImIT(p 2 ) can be 
approximated by ImII(r7i|j,) over the width of the resonance. If the resonance is broad, the 
full energy dependence of n(p 2 ) must be taken into account. 
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where = j&jr^Pa,b are the partial decay widths of the resonance in chan- 
nels (a,b). The common relation Imllfjo 2 ) = \fp 2 ~T(p 2 ) is not used to avoid 
confusion (the width of the resonance in the pole approach is related to the 
pole position, it is not given by the tree level result following from the optical 
theorem). The renormalized mass tur and the tree-level coupling constants 
of the resonance to the two channels g a ,b are the parameters to be fitted when 
confronting the scattering amplitude (fT2|) with data. Once these parameters 
are extracted from experimental data, the pole mass m p and pole width Y p 
of the resonance are obtained from the pole equation 

D(s p ) = s p -m 2 R + Ren+(s p ) - ReU + (m 2 R ) + ilmU + (s p ) = , (13) 

with s p = m 2 — im p T p and n + (s) = n(s-He). The pole equation (fT3*j) involves 
a complex function of a complex variable. If for real s, IT + (s) = R(s) + il(s), 
then for arbitrary complex s 

n+(s) = ReR(s) - ImJ(s) + i [ImR(s) + ReJ(s)] . (14) 

In order to find all the poles associated with a resonance coupled to channels 
(a, b) we have to look for the poles of Eq. ()13|) in the four different Riemann 
sheets defined by the complex channel momenta p atb . Following the con- 
ventional classification, the sheets are enumerated according to the signs of 
(Imp a ,Imp fe ): 



sheet I (+ 

sheet II (— 

sheet III (- 

sheet IV (+ 



(Imp a > 0, Ivnp h > 0) , 
{Impa < 0, lmp b > 0) , 

(15) 

(Imp a < 0, lmp b < 0) , 
(Imp a > 0, lmp b < 0) . 



We restrict ourselves to the case of two-body channels involving particles of 
the same mass. Thus, the thresholds for channels (a, b) are 2m a ,b and the 



momenta are defined as p a ,b{s) = J s — 4m 2 b /2 (m a < m b is assumed). Since 
the complete propagator D(s) is indeed an explicit function of the momenta 
p a; b, D(s) = D[s,p a (s),pb(s)}, a change in the sign of the imaginary part of 
the momentum — a change of sheet — is achieved with the replacement of 
Pa,b by —p a ,b in the propagator. Therefore, the poles are found solving the 
following four pole equations: 

D[s,p a (s),p b (s)] , D[s, -Pa(s),p b (s)} , 

(16) 

D[s,-p a (s),-p b (s)] , D[s,p a (s),-p b (s)\ . 
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For the case of interest, namely the / (980) scalar resonance coupled to a 
pair of pions and a pair of kaons 8 , the real and imaginary parts of the finite 
part of the 1PI two-point function Tl(s) are 



R(s) 




^[2-Arlog(^)]+%F[2 



^ lo s(i^)] e 



K 




+ 



%^[2-2^arctan(^-)]6 K , 



(17) 



where fa = y 1 — Amf/s for i = n, K, (3k = \J4m 2 K /s — 1, 0# = 0(5 — 4777^), 
and 0j<- = 0(4771^ — s). It is worth remarking that the step functions are 
not introduced by hand but result from the present calculation and play a 
crucial role in the determination of the pole structure. 

So far we have discussed the framework needed for the description of a 
resonance coupled to two channels. However, there are not current exper- 
iments that allow for a direct comparison of two-particle scattering ampli- 
tudes with experimental data. Therefore, in order to carry out the numerical 
analysis one has to rely on production processes such as J/if> — > ^(tttt, KK) 
and J/if> — > u;(7r7r, KK) or central production in proton-proton scattering 
pp — > p(7rn, KK)p. As stated in the Introduction, we will perform our 
analysis using only the former J /if) — > (f)(nir, KK) decays as a mechanism 
for producing pairs of pions and kaons. In this respect, we follow Ref. [Sj 
to relate the production amplitude F, also constrained by unitarity, to the 
scattering amplitudes T ab in Eq. (fT2|). The corresponding amplitudes for 
J /if) — > </>(7T7r, KK) are then written as 



where the real coupling functions a^^s) are parametrized as Oj(s) = 7^0 + 
7iiS and the 7« are obtained from the fit. Note that the theoretical expression 
for the mi phase shift and inelasticity are obtained from the T n7r scattering 
amplitude in Eq. (jHJ). The tttt data used in our fits are extracted from three 
different experimental analyses, two of them based on data from the 

reaction ir~p — > 7r + 7r~n [2H] and the third one fZT\ from 7c~p^ — > 7r + 7r~n 



8 In our analysis, we work in the isospin limit and therefore the mass difference between 
K° and K + is not taken into account for the KK threshold. The inclusion of this mass 
difference would deserve a more refined 3-channel analysis with eight Riemann sheets that 
is beyond the scope of the present work. 



F.,, = F(J/lf) -> 07T + 7T-) = 

F K = F(J/if) -> <f)K + K-) 





HH and n-p -> 7r°7r°n |3T]. 
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3 Numerical analysis 



Before proceeding with the numerical analysis we should keep in mind that 
our method is based on the pole approach and thus, in order to obtain a 
background independent fit to the data, we need to restrain ourselves to a 
neighborhood of the /o(980) resonance. For this reason we have chosen to 
work with experimental data on J/^ </>7r + 7r~ and J/ip — > <j)K + K~ decays 
[TBI E2] and on the tttt phase shift and inelasticity (251 1213 I2H EE] in the 
range 0.8 < y/s < 1.1 GeV. It is worth noticing that within the kinematical 
region between 0.8 GeV and the point where the rising of the tttt phase shift 
starts due to the appearance of the / (980), the contribution of the a scalar 
resonance to the phase shift in this region can be reasonably described in 
terms of an energy polynomial 9 (see below). 

For the tttt phase shift we have used three different sets of data. The 
first two sets differ mainly in the point lying just around 980 MeV and 
correspond to solutions B [23] (including this controversial point) and D 
j2B| (not including it) of Ref. [29J. The last set of data corresponds to the 
"down-flat" solution of Refs. (23 |2H] that seems to be the most preferable 
solution after a joint analysis of the S'-wave tt + tt~ and 7r°7r° data [30, 131]. 
In the following we will denote these three sets of tttt phase shift data as 
sets B, D and DF respectively. Part of the differences concerning the pole 
parameters of the / (980) resonance reported in the literature could be due 
to the use of different data. In this work, in order to quantify the influence 
of the tttt phase shift used, we have performed fits using data-sets B, D and 
DF. 

In the data fitting, a background term can be introduced for one or both 
channels and furthermore different energy dependences of the phase shifts 
can be considered. In our analysis we have included a background for the tttt 
and KK channels both with an energy dependence ranging from constant 
(5 = bo) to quadratic (5 = bo + bis + b 2 s 2 ). For the tttt case, the contribution 
of the a resonance to the background term is included in this way. For the 
KK case, only b\ and 62 are independent parameters since the background 
term must vanish by continuity below the kaon threshold. So then, the 
maximum number of parameters of our fits is 12: the renormalized mass 
rriR, the coupling constants gf 07T n and gt^KKi t ne parameters 6q, b\ and b\ for 
the background term for pions, bf- and b£ for the kaon background, and the 
constants 7^0,7^1 and 7ko ; 7xi parametrizing the F n and F K amplitudes. 

The results of the different fits performed show that i) the \ 2 improves 

9 Moreover, the choice of using experimental data only around 980 MeV avoids the 
need of describing the broad bump seen in the tttt phase shift around 600 MeV for which 
a simple polynomial parametrization is not adequate. 
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Fit 


set B 


set D 


set DF 


m\ (GeV 2 ) 


0.966 ±0.003 


0.982 ±0.001 


0.978 ±0.003 


g 2 fo7T JlQn (GeV 2 ) 


0.071 ±0.007 


0.065 ±0.005 


0.11 ±0.01 


91 K k/IQtt (GeV 2 ) 


0.25 ±0.03 


0.16 ±0.01 


0.31 ±0.04 


^0 


1.89 ±0.09 


1.47 ±0.03 


0.49 ±0.05 


67 (GeV" 2 ) 


-1.1 ±0.2 


1.47 ±0.09 


4.6 ±0.1 


6? (GeV" 4 ) 


0.8 ±0.1 


-0.50 ±0.06 


-2.83 ±0.09 


bf (GeV" 2 ) 


28.7 ±0.6 


15.5 ±0.3 


26.6 ±0.4 


b« (GeV" 4 ) 


-12.6 ±0.2 


-6.8 ±0.1 


-11.5 ±0.2 


7tt0 


5.1 ±0.4 


-0.2 ±0.9 


5.0 ± 1.0 


77rl (GeV" 2 ) 


-1.5 ±0.5 


8.2 ± 1.0 


1.7 ±1.0 




-27.2 ±0.5 


-39.5 ±0.7 


-21.7±0.5 


lKl (GeV- 2 ) 


29.7 ±0.5 


45.0 ±0.7 


26.2 ±0.6 


X 2 /d.o.f 


1.10 


0.87 


1.10 



Table 1: Values for the parameters obtained from a joint fit to the J/ip — > 
0(7T7r, KK) decays and the inr phase shift and inelasticity. Set B, D or DF 
refers to the set of nn phase shift data used in the fit. 

when a KK background is included, although changing its energy depen- 
dence makes no relevant difference; ii) the values obtained for the phys- 
ically relevant parameters (rriR, gf 07T n and gf Q Ki<) change only a few per- 
cent when different backgrounds are considered. The outcome of the fit 
for the case of quadratic backgrounds for pions and kaons is written in Ta- 
ble [0 and shown in Fig. The values for the renormalized mass m^, the 
coupling constants gf o7VK and g^KR-, the \ 2 P er degree of freedom and the 
other fitted parameters are presented for sets B, D and DF of tttc data. 
Concerning the values of the coupling constants obtained from the fit, we 
observe that the coupling of the / (980) to kaons is stronger than the cou- 
pling to pions: g 2 oK ^/ g 2 onn = 3.52,2.46,2.82 for sets B, D and DF respec- 
tively. Model dependent values for these coupling constants have been re- 
cently reported by the SND [33], CMD-2 [3Ij and KLOE [35] Collaborations. 
Their analyses, based on the study of the <fi — > 7r°7r°7 radiative process, give 
g 2 foK+K ./g 2 fon+n . = 4.6 ± 0.8 (SND), 3.61 ± 0.62 (CMD-2) and 4.00 ± 0.14 
(KLOE) 10 . Other analyses, based either on pp —>■ p(7C7r, KK)p central pro- 
duction [11J or on / (980) production in hadronic Z decay [3E], suggest the 

10 The coupling constants <?/ o7r7r (K /?) used in our analysis are related to the more common 

2 

coupling constants gf o7r +Tr-(K+K-) an d 9-k(k) used in experimental analyses by | 9, £™ = 

2 2 2 

g /o^+x- _ 4 2 J 1 g SgKK _ 9 f K+K- _ 2 

4^ - 3 9* 171 f and 2 ~ 4^ - 9K m fo . 
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same behaviour for the coupling constants and obtain gx/ 9-k — 2.1 ± 0.6 
and gx/g-K — 10 respectively. On the contrary, the analysis by the E791 
Collaboration |37j on the / (980) production in D s — > 3n decays gives 
gx/gn = 0.2 ± 0.6. The values we obtain for both coupling constants and 
for their ratio are smaller than model predictions [3*8*1 "4T)j and also than 
previous determinations [31]. 

Let us now proceed to determine the mass and width of the fo resonance 
within the pole approach. Fits to the data sets B, D and DF lead to values 
for the renormalized mass of uir = 983 ± 2 MeV, m# = 991 ± 1 MeV and 
rriR = 989 ± 2 MeV respectively. The pole parameters of the resonance are 
determined once the values of the renormalized mass and the ones regarding 
the coupling constants (see Table HJ) are included in the pole equation (fT3~j) . 
The numerical solution of the pole equation yields for data sets B, D and 
DF: 

= 987 ± 3 MeV , = 42 ± 9 MeV , 

m° = 999 ± 2 MeV , = 39 ± 8 MeV , (19) 

m® F = 1001 ± 6 MeV , T° F = 52 ± 16 MeV . 

This determination of the pole mass and width of the /o(980) resonance 
together with its couplings constants to the irn and KK channels constitute 
the main result of this work. Our results in Eq. (fT9*|) are in fair agreement 
with several values of the / (980) pole parameters appeared recently in the 
literature. The values m p = 994 MeV and T p = 28 MeV are obtained 
from an analysis of meson-meson interactions in a nonperturbative chiral 
approach [12! ■ Similar analyses give m p = 987 MeV and T p = 28 MeV [HI] or 
m p = 981.4 MeV and T p = 44.8 MeV [Sj. Other analysis based on the study 
of meson-meson interactions in different coupled channel unitarity models 
give m p = 1015 ± 15 MeV and T p = 86 ± 16 MeV gS] , m p = 991 ± 3 MeV 
and T p = 71 ± 14 MeV [H], rn p = 1008 MeV and T p = 54 MeV 0S|, 
m p = 993.2 ± 6.5 ± 6.9 MeV and r p ~ 100 MeV or m p = 1006 MeV and 
T p = 34 MeV @Bj. 

In addition to this result we have also analyzed the variation of the pole 
position as a function of the renormalized mass thr, i.e. the coupling con- 
stants g/ 07 r7r and g f KK are kept fixed to their values for set D in Table [TJ In 
Fig. [21 we show m p versus m p T p , which are related to the real and imagi- 
nary parts of the pole s p , for values of the renormalized mass in the range 
941 < rriR < 1027 MeV. Thus, each point on the plot corresponds to a solu- 
tion of the pole equation (|13[) — in terms of the values obtained for m p and 
T p — for a given value of m^. Only the physically relevant pole is shown in 
Fig. [21 complex conjugate poles or any other kind of poles are not included. 
In order to generate Fig. [21 we looked for solutions of the pole equation in 
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Dm2 pion mass spectrum GeV Dm2 kaon mass spectrum GeV 




0.9 0.95 1 1.05 1.1 1 1.025 1.05 1.075 1.1 

Mk3 pion mass spectrum GeV Mk3 kaon mass spectrum GeV 



Figure 1: Fit to tttt data on the phase shift and inelasticity (cross [25] , 
solid square [2S], solid triangle [2ZII2H]) and to data on J/^ ^ 0(7T7r, KK) 
production [THIISS]. Fits to the data sets B, D and DF are shown with dotted, 
dot-dashed and solid lines respectively. 
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0.9 0.925 0.95 0.975 1 1.025 1.05 



m 2 (GeV 2 ) 

Figure 2: Behavior of the pole position (m 2 versus m p T p ) as a function of 
the renormalized mass m# in the range 941 < uir < 1027 MeV. The set of 
points with larger values of |m p r p | appearing above 0.975 GeV 2 correspond 
to values of > 1020 MeV and are associated to poles in sheet III. 

each of the four Riemann sheets reaching the following conclusions: 

1. We did not find a pole neither in sheet I nor IV in the vicinity of 980 
MeV (we looked for poles in the range 960-1020 MeV). 

2. We find poles in sheet II in the range 941 < uir < 1027 MeV. From 
these solutions we see that the pole mass is always larger than the 
renormalized mass, i.e. m 1 ^ > ttir. When 4m 2 < m 2 R < Am 2 K , we 
obtain a pole width T p that has to be identified with the fo — > itn decay 
width. The width so found does not coincide with the decay width 
calculated from the tree level expression r^,-^ = gj o7T7T P n /167iy/s, the 
difference arising from the contribution of the IhlR(s) term in Eq. (JTlJl 
once s takes a complex value. 

3. We find poles in sheet III only for > 1020 MeV. In this case, the pole 
mass is always of the order of 20 MeV smaller than the renormalized 
mass, i.e. m p n < uir, and T p should be identified with the tree level 
width Tf ^ n7T + Tf _i K j{. Again, this decay width does not coincide 
with the pole width for the same reasons as before. 

From points 2 and 3 above we conclude that only one pole in sheet II will 
be necessary to describe the /o(980) resonance, since poles in sheet III only 



13 



appear for r/ijj > 1020 MeV while our fit to data always yields m# < 1 
GeV. Moreover, these points also indicate that the overlapping of poles is 
not possible. The idea of the overlapping of poles states that values of the 
renormalized mass niR close but below the kaon threshold {trr < 2m k) niay 
lead to values of the pole mass m p above the threshold (m p > Ivtik) and that 
values of m R close but above the threshold {m R > 2m k) may lead to values 
of m p below the threshold (m p < 2m K ). We stress that our conclusion is 
not general, it depends upon the values used for the coupling constants, in 
particular on the ratio g 2 foKR / g 2 fo7T7T gS]. 

4 Conclusions 

Using a generalized version of the Breit-Wigner parametrization based upon 
unitarity and a propagator obtained from effective field theory including loop 
contributions, we have performed a fit to experimental data on the tttt phase 
shift, the inelasticity and the J/ip — > 07T7T and J/ip — > (j>KK decays. The fit 
has been restricted to a neighborhood of 980 MeV (0.8 < y/s < 1.1 GeV) 
thus providing a process and background independent way of extracting the 
intrinsic properties (pole mass and width) of the /o(980) scalar resonance. 

The solution of the pole equation for the values of the parameters resulting 
from the fit allows us to conclude that the /o(980) is described in terms of a 
single pole in sheet II and yields the values mj = 987 ± 3 MeV and Tj = 
42 ±9 MeV, mf = 999 ±2 MeV and T% = 39 ±8 MeV, and mf o F = 1001 ±6 
MeV and Tf Q F = 52 ± 16 MeV for the data-sets B, D and DF respectively. 
We also analyzed the behaviour of the pole position as a function of the 
renormalized mass m R and found that a pole in sheet III only arises when 
mn > 1020 MeV while our fit to data yields always values m R < 1 GeV. 

Acknowledgements 

R. E. acknowledges A. Farilla, F. Nguyen, G. Venanzoni and N. Wu for 
valuable discussions, A. Bramon for a careful reading of the manuscript, 
L. Lesniak for providing us with the latest ittt data, and J. R. Pelaez for 
the Mathematica program for extracting the poles from the scattering am- 
plitudes. Work partly supported by the EU, EURODAPHNE (TMR-CT98- 
0169) and EURIDICE (HPRN-CT-2002-00311) networks. Work also partly 
supported by the Ministerio de Ciencia y Tecnologi'a and FEDER (FPA2002- 
00748), and CONACYT. 



14 



References 

[1] D. E. Groom et al. [Particle Data Group Collaboration], Eur. Phys. J. 
C 15 (2000) 1. 

[2] T. Bhattacharya and S. Willenbrock, Phys. Rev. D 47 (1993) 4022. 

[3] D. Morgan and M. R. Pennington, Phys. Rev. D 48 (1993) 1185. 

[4] B. S. Zou and D. V. Bugg, Phys. Rev. D 48 (1993) 3948. 

[5] D. Morgan and M. R. Pennington, Phys. Rev. D 48 (1993) 5422. 

[6] B. S. Zou and D. V. Bugg, Phys. Rev. D 50 (1994) 591. 

[7] A. Bernicha, G. Lopez Castro and J. Pestieau, Phys. Rev. D 50 (1994) 
4454. 

[8] M. Benayoun, S. Eidelman, K. Maltman, H. B. O'Connell, 
B. Shwartz and A. G. Williams, Eur. Phys. J. C 2 (1998) 269 
larXiv:hep-ph/9707509| . 

[9] M. Feuillat, J. L. Lucio M. and J. Pestieau, Phys. Lett. B 501 (2001) 
37 |arXiv:hep-ph/0010145 . 

[10] A. Bernich a, G. Lopez Castr o and J. Pestieau, Nucl. Phys. A 597 (1996) 
623 |arXiv:hep-ph/9508388| . 

[11] D. Barberis et al [WA1 02 Collaboration], Phys. Lett. B 462 (1999) 462 
|arXiv:hep-ex/9907055| . 

[12] D. Barberis et al [WA1 02 Collaboration], Phys. Lett. B 453 (1999) 325 
|arXiv:hep-ex/9903044| . 

[13] D. Barberis et al. [WA102 Collaboration], Phys. Lett. B 453 (1999) 316 
|arXiv:hep-ex/9903043| . 

[14] D. Barberis et al [WA1 02 Collaboration], Phys. Lett. B 453 (1999) 305 
|arXiv:hep-ex/9903042| . 

[15] R. Bellazzini et al. [GAMS Collaboration], Phys. Lett. B 467 (1999) 
296. 

[16] D. Aide et al. [GAMS Collaboration], Phys. Lett. B 397 (1997) 350. 
[17] J. E. Augustin et al. [DM2 Collaboration], Nucl. Phys. B 320 (1989) 1. 

15 



A. Falvard et al. [DM2 Collaboration], Phys. Rev. D 38 (1988) 2706. 



N. Wu, ||arXiv:he~ ex/0104050 

U. G. Meissner and J. A. Oiler, Nucl. Phys. A 679 (2001) 671 
|arXiv:hep-ph/0005253| . 

R. J. Eden, P. Landshoff, D. Olive and J. C Polkinghorne, "The analytic 
S'-matrix," New York (1966). 

R. G. Stuart, Phys. Lett. B 262 (1991) 113; 

A. M. Badalian, L. P. Kok, M. I. Polikarpov and Y. A. Simonov, Phys. 
Rept. 82 (1982) 31. 

M. E. Peskin and D. V. Schroeder, "An Introduction To Quantum Field 
Theory," Reading, USA: Addison- Wesley (1995) 842 p. 

B. Hyams et al, Nucl. Phys. B 64 (1973) 134 [AIP Conf. Proc. 13 (1973) 
206]. 

G. Grayer et al, AIP Conf. Proc. 13 (1973) 117. 

R. Kaminski, L. Lesniak and K. Rybicki. Z. Phys. C 74 (1997) 79 
|arXiv:hep-p'h /9606362] . 

R. Kaminski, L. Lesniak and K. Rybicki, Eur. Phys. J. directC 4 (2002) 
1 |arXiv:hep-ph/0109268| . 

G. Grayer et al, Nucl. Phys. B 75 (1974) 189. 

H. Becker et al. [CERN-Cracow-Munich Collaboration], Nucl. Phys. B 
151 (1979) 46. 

J. Gunter et al . [E852 Collaboration], Phys. Rev. D 64 (2001) 072003 
|arXiv:hep-ex/0001038| . 

U. Mallik, in Strong Interactions and Gauge Theories, Proceedings of 
the XXIst Rencontre de Moriond, Les Arcs, France, 1986, edited by 
J. Tran Thanh Van (Editions Frontieres, Gif-sur-Yvette, 1986), Vol. 2, 
p. 431; W. Lockman, in Hadron'89, Proceedings of the 3rd International 
Conference on Hadron Spectroscopy, Ajaccio, France, 1989, edited by 
F. Binon et al. (Editions Frontieres, Gif-sur-Yvette, 1989), p. 709. 

[33] M. N. Achasov et al, Phys. Lett. B 485 (2000) 349 
|arXiv:hep-ex/00050l"7| . 



16 



[34] R. R. Akhmetshin et al. [CMD-2 Collaboration], Phys. Lett. B 462 
(1999) 380 |arXiv:hep-ex/9907006| . 

[35] A. Aloisio et al. [KLOE Collaboration] , |arXiv:hep-ex/02040i3l 

[36] K. Ackerstaff et al. [OP AL Collaboration], Eur. Phys. J. C 4 (1998) 19 
|arXiv:hep-ex/98020T3| . 

[37] E. M. Aitala et al. [E7 91 Collaboration], Phys. Rev. Lett. 86 (2001) 765 
|arXiv:hep-ex/0007027| . 

[38] N. A. Tornqvist, Eur. Phys. J. C 11 (1999) 359 |arXiv:hep-ph/9905282| ; 
Erratum: Eur. Phys. J. C 13 (2000) 711. 

[39] M. Napsuciale, arXiv:hep-ph/9803396; M. Napsuciale and S. Rodriguez, 
Int. J. Mod. Phys. A 16 (2001) 3011. 

[40] E. P. Shabalin, Sov. J. Nucl. Phys. 42 (1985) 164 [Yad. Fiz. 42 (1985) 
260]; Erratum: Sov. J. Nucl. Phys. 46 (1987) 768. 

[41] R. Kaminski, L. Lesniak and B. Loiseau, Eur. Phys. J. C 9 (1999) 141 
|arXiv:hep-ph/9810386| . 

[42] J. A. Oiler, E. Oset and J. R. Pelae z, Phys. Rev. D 59 (19 99) 074001 
[Erratum-ibid. D 60 (1999) 099906] |arXiv:hep-ph/9804209| . 

[43] J. A. Oiler and E. Oset, Phys. Rev. D 60 (1999) 074023 
|arXiv:hep-ph/98093"37| . 

[44] J. A. Oiler and E. Oset, Nucl. Phys. A 620 (1997) 438 [Erratum-ibid. 
A 652 (1999) 407] |arXiv:hep-ph/9702314| . 

[45] V. V. Anisovich, Phys. Usp 41 (1998) 419 [Usp. Fiz. Nauk 168 (1998) 
481] |arXiv:hep-ph/9712504| . 

[46] M. P. Loch er, V. E. Markushin a nd H. Q. Zheng, Eur. Phys. J. C 4 
(1998) 317 |arXiv:hep-ph/9705230| . 

[47] S. Ishida, M. Ishida, H. Takahashi, T. Ishida, K. Takamatsu and 
T. Tsuru, Prog. Theor. Phys. 95 (1996) 745 |arXiv:hep-ph/9610325| . 

[48] N. A. Tornqvist and M. Roos, Phys. Rev. Lett. 76 (1996) 1575 
[arXiv:hep-ph/95112To| . 

[49] R. J. Eden and J. R. Taylor, Phys. Rev. Lett. 11 (1963) 516. 



17 



